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Abstract 



> 

It is shown that the main contribution to the rotational curve of a spiral 
! galaxy may be due essentially to the interaction, in the general relativistic 

spacetime, of the galactic matter with a very light long range scalar field which 
respects the weak equivalence principle. The comparison of the theoretical 
results with 23 spiral galaxy rotation curves shows a good agreement between 
q ' our proposal and observations. 

cr 

99' 1 Introduction 



According to newtonian dynamics and Newton's inverse-square force law of gravita- 
tion, the circular velocity around an isolated body should decrease with the distance, 
r, to the centre like However, observations show that the rotation curves of many 
spiral galaxies flatten at large distances. The current interpretation consists to in- 
voke the existence of dark matter though its nature and distribution in space are 
not yet clearly defined [1]. 

Alternative proposals, which provide explanations for observed galactic rotation 
curves without the need of dark matter, have been developed by several authors. 
In that line of thinking, the Modified Newtonian Dynamics (MOND) proposed by 
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Milgrom [2] introduces a " universal" constant ao of the dimension of an acceleration 
such that newtonian dynamics breaks down at acceleration much smaller than ao (or 
equivalently Newton's law of gravitation fails when the magnitude of the potential 
gradient is much smaller than ao). This acceleration constant turns out to be of 
the order of cH suggesting a cosmological link [3]. However, actually the equality 
a = cH has no theoretical basis in the framework of MOND and it turns out 
that at the cosmological level MOND fails to satisfy the cosmological principle [4]. 
Besides, while having some successes at the scale of galaxies [5, 6, 7] and clusters 
of galaxies [8, 9], the MOND fails in the laboratory 1 [10]. Moreover, at the scale of 
clusters of galaxies, Gerbal et al. have found that a certain amount of dark matter 
is required even with MOND [11, 12]. 

In the framework of the conformal Weyl gravity, Mannheim and Kazanas ob- 
tained a complete, exact exterior solution for a static, spherically symmetric source [13] 
In addition to the exterior Schwarzschild solution, their solution contains an extra 
gravitational potential term which grows linearly with the distance r to the cen- 
tre. A cosmological link is also made by these authors who noted intringuingly that 
their parameter 7 is roughly the value of the inverse of the Hubble length, -^k Their 
solution is able to interpolate between those of Robertson- Walker and Schwarzchild 
in a continuous and smooth manner by exploiting the conformal structure of the 
conformal Weyl gravity. This solution has been applied successfully by Mannheim 
to the galactic rotation curves of spiral galaxies [14]. However, by matching the 
exterior solution to an interior one that satisfies the weak energy condition and a 
regularity condition at the centre, Perlick and Xu [15] have shown that this leads 
to contradiction of Mannheim and Kazanas's suggestion. They conclude that the 

conformal Weyl gravity is not able to give a viable model of the solar system. Be- 

1 This conclusion is rejected by M. Milgrom, since the experiment has been performed on earth 
where the gravitational pulling is much greater than ao- 
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sides, it seems that the cosmological models derived from the conformal Weyl gravity 
fails to fulfill simultaneously the observational constraints on present cosmological 
parameters and on primordial light element abundances [16]. 

Other authors have suggested that Newton's law of gravity, which describes so 
well the motion of moons and planets in the solar system, may break down over 
distances comparable to the size of a galaxy. So, alternative force laws have been 
proposed [17, 18, 19]. However, one should bear in mind that newtonian mechanics 
is the only appropriate approximation of general relativity in the weak field and low 
velocity limit (see the criticism of D. Lindley, [31]). 

In this paper, we propose a new dark matter scheme in the framework of Ein- 
stein's general relativity and its weak field and low velocity limit, newtonian gravity. 
Our new dark matter candidate is a long range neutral massive scalar field directly 
coupled to matter, unlike the Brans-Dicke theory in which the scalar field does not 
exert any direct influence on matter (its only role is that of participant in the field 
equations that determine the geometry of spacetime) [20]. In addition, one requires 
that the scalar field under consideration, 0, respects the weak equivalence principle. 
Let us notice that, if one takes into account the gravitational coupling of 0, with the 
ordinary matter, then the g a p s will depend not only on the x^ s but also on 0, that 
is 9a/3 = 9a/3(x^, 0). The plan of this paper is as follows. In section 2, we establish 
the equation of motion of a test body in the presence of the light scalar field, 0. In 
section 3, we determine the potential that couples the ordinary matter to the 
scalar field. In section 4, the results of both sections 2 and 3 are combined and 
applied to the rotational curves of spiral galaxies. Finally, in section 5, twenty three 
rotation curves are confronted with the results of the previous sections. 
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2 The equation of motion in the presence of a 
long range scalar field 

In the absence of any external field, the equation of motion of a test body writes in 
general relativity 

{u a V a )u» = (1) 

where u a = ^f- denotes the velocity four-vector of the test body and ds 2 = 
g a fjdx a dx 13 is the metric and V a is the covariant derivatives (a, f3 = 0,1,2,3). 
In the same manner as in the presence of an electromagnetic field, in the presence 
of a scalar field, <p, a force term enters in the right-hand side of equation (1). This 
force term would be a priori of the simple form — jjz^k involved by the lagrangian 
L = —\g a pu a vP — ^r, where the potential V<p depends on (V^ = if <p = 0) and is 
assumed to be the same for all matter at any point of spacetime to ensure obedience 
to the weak equivalence principle. Thus one would write 

(«^K = -^. (2) 

However, equation (2) is not satisfactory, since the unitarity condition u^u^ = 1 
implies 

u»(u a V a ) UlI = 0. (3) 

In order to satisfy relation (3), there should be at least the additional term js^-u^ 
to the right-hand side of equation (2). Finally, the correct equation of motion writes 



1 dVj ]_dV„ 

c 2 dx^ c 2 

or making appear explicitly the Christoffel symbols, 



«W^ ? ^ + ^v (4) 
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One may derive equation (4, 5) from the lagrangian 



L = -_(^ Q/JU V + 1). (6) 

2.1 Small speeds and weak fields approximation 

For small speeds and weak gravitational fields, with respect to a galilean referen- 
tial, it comes « -\d a g m . Now, g a p(x^,(/)) « ^(a;'*, 0) + (%f) 0=Q since we 
are interested by the weak 0-field approximation. One may also write g a p{x^,(f)) ~ 
(7 a/ g(x M ,0) + (7^) V^, for Vtf, is proportional to (f) at the first order. More- 
over, for small speeds and weak gravitational fields g 00 (x^,0) « 1 + 2^- whence 
S'ooC^^, 0) ~ 1 + 2 VN+ ^ 1 J~^ v ' t ' , where Vat denotes the newtonian gravitational poten- 
tial and we have set (f^) y = ^(1 + /) (h will be found in a next study on the 
mass distribution that / is positive and typically in the range 1CT 9 — 1CT 3 in spiral 
galaxies). Consequently equation (5) reduces at the first order to 

^_v (M+ if^ (7) 

where ^ = ^ + ^f + ^f + ^^in spherical coordinate (r, 0, <p). Particularly, 
for planar orbits in plane polar coordinates (r, 6), one gets from the above equation 

dv L _v[ d(V N + fV£ 1 dV^ 

dt r~ Or + C 2 dt Vr [ ' 

1 d{v e r) _18(Vn±J%1 IdV^ (q) 
r dt r 86 c 2 dt V ° [ } 

where v r = ^ is the radial velocity, vq = the tangential velocity. 

2.2 Static spherically symmetric fields 

For static spherically symmetric fields, equations (8) and (9) become 
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_v], = _ G M£ _ ev, 

at r r 2 or 

1 d(v e r) 1 dVtp 

since v r <C c, ^ = = and = G^p-, where G is the gravitational constant 
and m(r) denotes the mass up to radius r. On integrating through equation (11), 
one gets 

exp m 

v e = J^^- (12) 
r 

where J is a constant which would represent the angular momentum per unit mass 
if the scalar field (p were not present. 

3 Derivation of the potential 

In as much as we are interested by the weak field approximation only, the potential 
is simply a linear function of the scalar field, </>. Now, the equation of a scalar 
field, 0, writes in the weak field approximation 

^0+(^) 2 = -y0. (13) 

The above equation, as it is well known, follows from the lagrangian density of a 
real scalar field L = \d^4>d^4> - U ((/>), where U(<f>) = \ {^f f (f> 2 + J V(pd(P denotes 
the potential energy of the scalar field <fi and is its mass. Hereafter, we assume 
a phenomenological effective potential energy of the form 

U(<P)= 1 -(^fn 2 + qr (14) 

where p and q are real constants. Of course, whenever p is an integer, p > 3. In 
the following we determine only the exponent p, in as much as the potential U((f>) 
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is assumed negligible with respect to the ordinary matter mass-energy. Besides, 
since we have considered a long range scalar field 4>, we may drop the mass term in 
equation (13) and consider the massless field equation 



dPd^ = -V(f). (15) 
Then, the static spherically symmetric solution = <f>(r) satisfies the following 



equation 

d 2 



+ zZ = V* (16) 



dr 2 r dr 

and the potential energy reduces to the power law 

U{<t>) = qcjf. (17) 
So, we look for a solution of the form 

0ocr fc . (18) 
Then replacing by r k in equation (16) above it comes 

y^Jt±R (19) 
and eliminating r in relation (19) above by replacing r by (p 1 ^ gives 

V oc k{k+ \)<P' 2 / k (20) 

UocJvh^y^j^ 1 ' 1 ^- ( 21 ) 

Therefore, relation (21) together with relation (17) gives 
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k=- T —r. (22) 

2 

It is natural to expect the exponent p to be a positive integer. Particularly, we may 
consider the special case where this exponent is even. Therefore, as p > 3, it comes 
p = 2(n + 1) where the integer n runs from 1 to oo. This yields oc r~ x l n and then 
the following potential 

V 4> = Kr' 1/n (23) 

where K is a real constant the magnitude of which depends on the strength of 
the interaction of the scalar field <f> with the ordinary matter in a given region of 
spacetime. 

4 Application to the rotational curves of spiral 
galaxies 

If a long range scalar field such that the one considered in this paper does exist, then 
as a dark matter candidate it may influence significantly the dynamics in a spiral 
galaxy and thus modify the shape of its rotation curve. Mostly, this may be done 
without the need of a great amount of dark matter of this kind. Indeed, setting 
™2 } rfc( r ) = fl^i, equation (9) takes the form 

_ v$ = c m(r) + mj ark (r) 
dt r r 2 1 ' 

For K < 0, the derivative ^ = — ^J™ n is positive and thence rn^fj rk (r) = /^H?" 1 ~ 1 ' /n 
mimics a dark matter mass profile which may be important though the energy cor- 
responding to the potential of the scalar field, U(<f>), is in fact rather negligible with 
respect to the real matter mass m{r). It is worth noticing that our interpretation 
of mlfj rk (r) as a hidden mass term involves that the integer n should run from 2 
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to oo (instead of 1 to oo) because one will require the positivity of the derivative 
^ in addition to the positivity of m^ rk (r) that is (r) = ^/^r" 1 /™ > 
everywhere within a galaxy (except at the boundary). In this section, we wish to 
prove the validity of relation (12) on the basis of experimental data available on the 
rotation curves of spiral galaxies. To test relation (12), it is better to rewrite it as 
follows 

In (rv ) = In J + ^ (25) 

and plot In (rv e ) versus r for each rotation curve, where r is the distance to the 
centre. More precisely we use the least-squares fit and search the coefficients a and 
b such that In (rvg) = ar~ l l n + b, where a = ^ and b = In J. It turns out that 
the potential is indeed well represented by a function of the form Kr~ 1 / n , the 
exponent n being the integer that leads in absolute value to the highest correlation 
coefficient, R, for a given rotation curve. Consequently, the coefficient a should 
be negative in as much as the interpretation of m^ rk {r) as a missing mass term 
holds. Table 1 below summarizes the numerical results obtained for twenty three 
spiral galaxies including three dwarf galaxies (DDO 170, NGC 3109 and the dwarf 
"regular" UGC 2259) and three giant low surface brightness disk galaxies (Malin 1, 
NGC 7589 and F 568-6). Clearly, the coefficient a is always negative whereas b is 
always positive. 

5 Fits to individual rotation curves 

The following figures show that the rotation curves of spiral galaxies (covering a wide 
range of mass and size) can be reproduced essentially by a potential, V^, resulting 
from the interaction of a long range scalar field with the usual galactic matter in as 
much as this interaction, like gravity, respects the weak equivalence principle. It is 
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found that fairly good fits can be obtained in the first approximation without the 
need to take into account the photometric properties of the galaxies. 



name [reference] 


n 


R 2 


a 


b 


DDO 170 [21] 


3 


0.999655088 


-7.23 


9.8 


DDO 170* [21] 


3 


0.99970403 


-7.1 


9.7 


NGC 3109 [22] 


10 


0.998736961 


-30.2 


26.29 


NGC 3109* [22] 


10 


0.998784224 


-30.17 


26.3 


UGC 2259 [23] 


7 


0.999719756 


-15.14 


16.89 


Malin 1 [24] 


3 


0.999950275 


-11.67 


12.33 


NGC 7589 [24] 


9 


0.99981229 


-16.11 


19.72 


F 568-6 [24] 


9 


0.999709336 


-17.67 


20.97 


NGC 4419 [25] 


6 


0.999743147 


-8.645 


13.486 


NGC 1035 [25] 


5 


0.99941517 


-8.4 


12.51 


NGC 1325 [25] 


13 


0.999646177 


-18.94 


23.3 


NGC 4062 [26] 


18 


0.999703876 


-24.49 


29 


NGC 2742 [25] 


15 


0.999726171 


-20.7 


25.2 


NGC 3067 [25] 


5 


0.99959742 


-7.58 


12.04 


NGC 247 [27] 


11 


0.998423849 


-25.3 


25.16 


NGC 3198 [28] 


5 


0.999638676 


-6.716 


11.52 


UGC 12810 [25] 


11 


0.999759876 


-16.57 


20.99 


NGC 4051 [29] 


8 


0.998352026 


-12.56 


16.61 


UGC 3691 [26] 


7 


0.999209314 


-12.57 


16.14 


NGC 3593 [25] 


2 


0.999652662 


-2.235 


6.91 


NGC 2639 [25] 


2 


0.998720197 


-5.77 


9.92 


NGC 4378 [25] 


6 


0.999521022 


-7.79 


13.34 


NGC 4448 [25] 


4 


0.998699541 


-6.306 


11.073 


NGC 7606 [25] 


3 


0.998640104 


-7.4 


11.3 


NGC 4402 [30] 


13 


0.99992101 


-27.85 


29.22 



Tab. 1. The exponent n, the square of the correlation coefficient, R 2 , and the 
parameters a and b of the least-squares fit to the rotation curves of 23 spirals are 
tabulated. The asterisk points out that the parameters presented are those of the 
final rotation curve obtained after a correction for asymmetric drift (caused by a 
pressure parameter) has been applied to the individual velocity points. As expected, 
the correlation coefficient increases when the correction for asymmetric drift is ap- 
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plied. 



6 Concluding remarks 

The work which has been presented in this paper suggests the existence of a neutral 
scalar particle which may contribute significantly to the rotation curves of spiral 
galaxies. The results involved by this proposal are found to be in good agreement 
with the available observational data. However, we have assumed that the referential 
attached, to the dynamical centre of a given galaxy is galilean. But, spiral galaxies 
are disks in differential rotation. So, a small correction may be necessary to fit with 
more accuracy some peculiar rotation curves. Besides, we have fitted the rotation 
curves of 23 spiral galaxies without the need of photometric data since in our study 
the main contribution is brougt by the scalar field. Nevertheless, this is not necessary 
problematical because it just brings one more support to the fact that there exists a 
strong correlation between the spatial distribution of the luminous matter and the 
spatial distribution of the dark matter in spiral galaxies. 
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